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We investigate the structure of magnetic field amplified by turbulent velocity 
fluctuations, in the framework of the kinematic Kazantsev-Kraichnan model. We 
consider Kolmogorov distribution of velocity fluctuations, and assume that both 
Reynolds number and magnetic Reynolds number are very large. We present 
the full numerical solution of the model for the spectra and the growth rates of 
magnetic fluctuations. We consider astrophysically relevant limits of large and 
small magnetic Prandtl numbers, and address both helical and nonhelical cases. 
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1. Introduction 

Magnetic fields are found everywhere in the universe, in planets and stars, in galaxies 
and galaxy clusters. Magnetic fields in astrophysical systems are usually generated in a 
broad interval of scales, ranging from small resistive scales to large scales exceeding the 
correlation length of plasma motions. One of the most important, challenging and still open 
questions in astrophysics is how cosmic magnetic fields have been generated and what is 
their structure. The prevailing theory for the origin of magnetic fields is dynamo action, 
which is stretching of magnetic field lines by r andom motion of highly conducting p 



asmas 



or fluids in whic: 


1 these lines are frozen (e.g.. 


Vainshtein & Zeldovich 


1972 


Parker 


1979; 


Lvnden-Bell 


1994 


Zwei 


3el & ] 


ieiled 


1997; 


Brandenburg & Subramanian 


2005; 


Kulsrud 


2005; 


Schekochihin & Cowlev 


2006; 


Kulsrud & Zweibel 


2008). Small-scale ma 


gnetic fields, that is 



fields at the scales smaller than the velocity field scales, are generally expected under this 
mechanism, while large-scale magnetic fields, correlated at scales larger than the correlation 
scale of a velocity field, can be generated if some additional conditions are met, such as 
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the condition that the velocity field v(x, t ) is not mirror symmetric, say, possesses nonzer o 



kinetic hehcity H = J v ■ {V x v) d^x ^ fISteenbeck. Krause fc Radlerlll966l : lMoffattlll978f ). 



Assume that the velocity fiuctuations have the typical correlation scale Iq, and the typ- 
ical rms value Vq, and the plasma has kinematic viscosity u and magnetic diffusivity rj (the 
latter is proportional to electrical resistivity). The range of scales available for velocity and 
magnetic fiuctuations can be characterized by the Reynolds number Re ~ loVf)/v and the 
magnetic Reynolds number Rm ~ IoVq/ti, respectively. In astrophysical applications both 
numbers are very large (for example, in a protogalaxy, where dynamo action is believed to 
operate. Re and Rm reach values ~ 10^ and ~ 10^^, respectively). Their ratio, the magnetic 
Prandtl number Pm = Rm/Re can be either large or small. For example, in galaxies and 
galaxy clusters Pm ^ 1, while in planets and stellar interiors Pm <^ 1. Because of the 
vast range of scales available for magnetic and velocity fiuctuations, and generally strongly 
disparate magnetic and velocity dissipation scales, present-day direct numerical simulations 
cannot directly address astrophysical magnetohydrodynamic (MHD) regimes. Indeed, max- 
imal Reynolds and magnetic Reynolds numbers accessible with numerical simulations are 
hopelessly small, of the order of 10^ — 10^, in which case the typical magnetic Prandtl num- 
bers are not significantly different from Pm ~ 1. A physical picture of magnetic dynamo 
action and an effective analytic framework for investigating astrophysical dynamo action are 
therefore in demand. 

The first step in understanding dynamo action is to understand the initial, kinematic 
stage of magnetic field amplification. In this regime, the magnetic field is weak and does not 
affect the velocity fluctuations. Magnetic fleld evolution is therefore fully described by the in- 
duction equation in which the velocity fleld is prescribed independently of the magnetic fleld. 
An effective framework in this case is provided by the so-called Kazantsev-Kraichnan model, 
where the velocity fleld is assumed to be a random Gaussian short-time-correlated fleld. This 
formal simpliflcation allows for analytic solutions of the model while capturing the essential 
physics of the phenomenon. We should note however that even with this simpliflcation the 
model is nontrivial and its general solution is not known. Only certain special cases have been 
solved so far, which reveal a good agreement with numerical simulations in the parameter 



range accessible to numerics (e.g ., iMaron &: Cowleyl l200ll : iHaugen. Brandenburg fc Dobler 



2004J : iBoldvrev fc Cattaneol 120041 ) 



The Kazantsev-Kraichnan dynamo model allows one to answer the fundamental ques- 
tions concerning the possibility of turbulent dynamo action for given Reynolds and mag- 
netic Reynolds numbers, the spectrum of growing magnetic fluctuations, the conditions for 
large-scale magnetic fleld ampliflcation, etc. In many instances, the results obtained in high- 
resolution numerical simulations were predicted by the model well before such simulations 
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become available. We therefore believe that the model can provide a valuable insight into 
astrophysical dynamo regimes that can hardly be accessed through direct numerical simula- 
tions in foreseeable future. 

Our preliminary results aimed at the full numerical characterizati on of the kinematic dy- 



namo action in the Kazantsev-Kraichnan framework were presented in lMalyshkin fc Boldyrev 



(120071 ). In particular, we investigated the growth rates of magnetic fluctuations in the veloc- 
ity field with the Kolmogorov energy spectrum. In the present paper, we explore the spatial 
structure of the growing magnetic eigenmodes. We assume the Kolmogorov spectrum of 
velocity fluctuations and address extremely large Reynolds numbers, up to Re ~ 10^. This 
allows us to study astrophysically relevant limits of very large and very small magnetic 
Prandtl numbers. For completeness, we perform our analysis for the cases of both large 
and small kinetic helicity. In the helical case we also discuss the relevance of the conven- 
tional a-model for the description of large-scale dynamo action, in the case when small-scale 
magnetic fluctuations are amplified as well. 

In the next section, we describe the Kazantsev-Kraichnan model of kinematic dynamo 
action. In Section 3, we present our results, and in Section 4 we give our conclusions. 



2. Kazantsev-Kraichnan model 



Kinematic dynamo action is described by the induction equation for the magnetic field: 

= V X (v X B) + r^V^B, (1) 

where v(x, t) is the velocity field, B(x, t) is the magnetic field, and rj is the magnetic diffu- 
sivity. In this equation the vel ocity field is presc ribed independently of the magnetic field. 
Following iKazantsevI (119681 ) and iKraichnanI (119681 ). we assume that the velocity field is statis- 
tically homogeneous and isotropic and has a Gaussian distribution with zero mean, (v) = 0, 
and the following covariance tensor 



{v\^, ty{^', t')) = (x - x')5(t - f), 
where is an isotropic tensor of turbulent diffusivity. 



(2) 



x^ 



+ K,L- 



XX' 



2 + ge'^'^x''. 



X 



(3) 



Here, functions kl(x) and g{x) describe kinetic energy and helicity, x = |x|, brackets () 
denote ensemble average, e*-'^ is the unit antisymmetric pseudotensor and summation over 
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repeated indices is assumed. The first two terms on the right-hand side of Equation ([3]) 
represent the mirror- symmetric, nonhehcal part, while the last term describes the helical 
part of the velocity fiuctuations. For an incompressible velocity field (the only case we 
consider here), we have kn{x) = kl{x) + xk'i^{x) /2, where the prime denotes derivative with 
respect to x = |x|. Therefore, to describe the velocity field, we specify only two independent 
functions, kl{x) and g{x). The Fourier transformation of Equation ([3]) is 



K'^{k) = F{k) [5'^ - 



(4) 



Functions F{k) and G{k) can be obtained from functions njjx) aiid q{x \ and vice versa, by 
using the three-dimensional Fourier transforms (IMonin fc Yaglonj|l97ll ). 



The correlator of homogeneous and isotropic magnetic field can similarly be expressed 



as 



(i?*(x,t)i?^(0,t)) = M^U^^- 



x'^x^ 



x^ 



+ Ml 



XX- 



+ Ke'^^xK 



x^ 



(5) 



where the field solenoidality constraint div B = implies Miy{x, t) = Ml{x, t) + (x/2)M^(x, t). 
To fully describe the magnetic field correlator, we need to find only two functions, ML{x,t) 
and K{x,t), corresponding to magnetic energy and magnetic helicity. The Fourier trans- 
formed version of Equation ([5]) is 



{B\k,t)B*^ik,t)) = Fnik,t){5'^ 



k'y 



(6) 



where Fsi^k, t) is the magnetic energy spectral function, (|B(k, = 2FB{k, t), and Hsik, t) 
is the spectral function of the electric current helicity, (S**(k, t) ie'^^'-k^ B^(k,t)) = HB{k,t). 
The problem is then to find the correlation function (Equation ([5])) of the magnetic field, or, 
alternatively, its Fourier version (Equation ([6])). 

Suppose that the velocity field (Equations ([2]) and ([3])) is given, i.e., kinetic energy kl{x) 
and kinetic helicity g{x) are given. In this case, to find the properties of the growing 
magnetic field driven by helical dynamo action, one needs to solve two coupled partial 
differential equations for functions Mi(x, t) and K(x, t) related to magnetic ene r gy an d 
magnetic helicity. Such equations were first derived by IVainshtein fc Kichatinovl (119861 ). 
Due to their complexity, there have been few theoretica l resu l ts obtained for the helical 
dynamo (e.g.. iKulsrud fc Anderson! Il992l: iKim fc Hughes! 119971: iBlackman fc Field 1 12002| ). 



Brandenburg fc Subramanian 



20^ and references 



For direct numerical simulations see 
therein. Recently, it has bee n established in iBoldyrev. Cattaneo &: Rosner! (12005! ) that 
Vainshtein &: Kichatinov! (11986! ) equations also possess a self-adjoint structure, which is simi- 
lar to a two-component quantum mechanical "spinor" form with imaginary time. These two 
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coupled self-adjoint differential equations are linear and homogeneous, and they describe 
the growth of the magnetic field in the Kazantsev-Kraichnan model with nonzero kinetic 
helicity. We solve them numerically by the fourth-order Runge-Kutta integration method 
and by matching the n umerical solution to the analy tical asymptotic solutions at x — > and 
x — > cx), for details see iMalyshkin fc BoldyrevI (l2007l ). 

We are interested in fast exponentially growing eigenmodes of the magnetic field am- 
plified by helical kinematic dynamo. Therefore, both magnetic correlator functions ML{x,t) 
and K{x, t) are assumed to be propor tional to exp(At), where A is the gro wth rates of the 
field eigenmodes. It is important that iBoldyrev. Cattaneo &: Rosnerl (120051 ) helical dynamo 
equations, which we solve, are self-adjoint because this guarantees that all growth rates 
A are real. It turns out that in analogy with quantum mechanics, there are two types of 
magnetic field eigenmodes: bound (spatially localized) and unbound (spatially nonlocal- 
ized). First, for growth rates A > Aq = (7^(0)/[/tL(0) + 2rj\ the eigenf unctions are bound and 
correspond to "particles" trapped by the potential provided by velocity fluctuations. The 
bound eigenmodes have discrete growth rates, i.e., A = A„ > Aq where n = 1,2,3.... The 
bound eigenf unctions decline exponentially to zero as x ^ oo. Second, for A < Aq the eigen- 
functions are unbound and correspond to "travehng particles". The unbound eigenmodes 
have continuous eigenvalues of their growth rates, < A < Aq. The unbound eigenfunctions 
asymptotically become a mixture of cosine and sine standing waves as x — > cxd. Eigenvalue Aq 
corresponds to the fastest growing unbound eigenmode. The structure and other properties 
of the magnetic field amplified by dynamo action are fully determined by all growing eigen- 
modes of the magnetic field. In particular, the magnetic energy spectral function FB{k,t) is 
the sum over all energy spectral eigenfunctions, 

"max ^ Ao 

(|B(k,t)p)/2 = FB(A;,t) = Vc„FB,„(fc,t)+ / c(A)Fb,a(A;, t) rfA. (7) 

n=i Jo 

In this equation Fs^nik, t) and FB,x{k, t) are the energy spectral eigenfunctions for the bound 
and unbound eigenmodes, respectively; coefficients c„ and c(A) depend on the seed magnetic 
field at the initial moment when dynamo started to operateJll Similarly to Equation ([7]), 
the electric current helicity spectral function Hsik^t) is the sum over all helicity spectral 
eigenfunctions with the same coefficients as in Equation ([7]). To find the correlation function 
of the magnetic field, it is sufficient to find all growth rates and the corresponding spectral 
eigenfunctions. 



^ Note that while the energy spectrum FB{k,t) = (|B(k, t)p)/2 stays always positive, the individual 
energy spectral eigenfunctions can be negative. The sharp declines in individual eigenfunctions FB^k) and 
Hsik) plotted in Figures [2] and [3] are due to the logarithmic representation of the eigenfunctions in these 
figures, which has difficulty to suit with values Fsik) — and Hsik) — 0. 
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Fig. 1. — Growth rates A„ of the bound magnetic eigenmodes, and Aq of the fastest growing 
unbound eigenmode. Plots (A) and (B) are for /i = 1 and 0.1, respectively, while fcmax = 3000 
(Reynolds number Re ^ 1 and Prandtl number Pm ^ 1). Plots (C) and (D) are for h = 1 
and 0.1, respectively, while fcmax = 3 x 10^ (Re ^ 1 and Pm <^ 1). All plots are on the 
logarithmic scale. 



3. Results 

To study a case relevant to astrophysical systems, we consider velocity correlation tensor 
(Equation (jl])) with the Kolmogorov power velocity spectrum and large Reynolds number. 
Re 1. In the Kolmogorov turbulence, the t urbulerit diffu sivity, given by Equation ([3]), 



scales as vil ~ l'^^^, where / = |x — x'| (see, e.g., lFrischlll995l ). As a result, the Kolmogorov 
scaling implies kl(x) ~ /t(0)(x//o)'^^^ ~ Volo{x/lo)*^^. Without loss of generality we take 
/o ~ 1, fo ~ 1, and therefore 

F{k) = k~^^/^, 

G{k) = -hk-^Fik), 2<A;<A;^ax. (8) 

Here, the lower cutoff wavenumber is fcmin = 2, the upper cutoff wavenumber fcmax ~ 
2[kl{0)/i']^^^ ~ 4i/~^/^ is determined by the plasma kinematic viscosity u, and the helic- 
ity parameter h must satisfy the realizability condition —l<h<l] the velocity field is 
maximally helical when \h\ = l|^ 

We choose magnetic diffusivity to he rj = 10"^, corresponding to magnetic Reynolds 
number Rm ~ 10^. By changing the value of the turbulence cutoff wavenumber fcmax in 



^ Given kl (x) , function q(x) cannot be chosen arbitrarily, its Fourier image must satisfy the reahzabihty 
condition |G(fc)| < F{k)/k ( Moffatt IQTSh . This results in the condition — 1 < /i < 1. Analogously, given 
ML{x,t), function K{x,t) is restricted by condition \HB{k,t)\ < 2kFB(k,t). 
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Equation ([8]), we vary the Reynolds number Re ~ 1/u ~ and the magnetic 

Prandtl number Pm = Rm/Re v/rj. We study two cases for the Prandtl number: a case 
when it is large and a case when it is small. In the first case, Pm ~ 150 ^ 1, which is 
achieved by choosing fcmax = 3000 (Re ~ 6800). In the second case, Pm ~ 6.7 x 10~^ <^ 1, 
which corresponds to our choice fcmax = 3 x 10^ (Re ~ 1.5 x 10^). These two cases of large and 
small Prandtl numbers are considered in combination with two cases for the kinetic helicity: 
first, a case when h = I'm Equation ([8]) and the velocity field is maximally helical and second, 
a case when /i = 0.1 and the kinetic helicity is small. Thus, in total we consider four cases 
for our choice of the Prandtl number Pm and the kinetic helicity parameter h. The resulting 
growth rates A„ of the bound (localized) eigenmodes and Aq of the fastest growing unbound 
(nonlocalized) eigenmode are shown on the logarithmic-scale plots in Figure [H The growth 
rates are measured in the units of large-scale eddy turnover rate ~ Vq/Iq {vq/Iq ~ 1 here). 
The logarithmic-scale plots of the absolute values of magnetic energy spectral eigenfunctions 
are given in Figure [2l The logarithmic-scale plots of the absolute values of electric current 
helicity spectral eigenfunctions are given in Figure [31 

In the case Pm ~ 150 and h = 1 there exist 10 growing bound eigenmodes of magnetic 
field, whose growth rates are shown on the plot (A) in Figure [H Among these we select 
four bound modes Ai ~ 35.41, A2 ^ 42.74, A7 ~ 213.4, and Aio ^ 731.1, and we plot their 
magnetic energy and current helicity spectral eigenfunctions by the dotted, dash-dotted, 
dashed, and smooth solid lines, respectively, on the left-upper plots in Figures [2] and [3l The 
spectral eigenfunctions of the fastest unbound eigenmode Aq — 33.23 are shown by the red 
jagged spiky lines on these left-upper plots. In the case Pm ~ 150 and h = 0.1 there are 
just four bound magnetic field eigenmodes, which are shown on the plot (B) in Figure [1] 
These eigenmodes are Ai ~ 0.3336, A2 ^ 26.03, A3 ^ 190.1, and A4 ~ 654.2, and their 
energy and helicity spectral eigenfunctions are shown by the dotted, dash-dotted, dashed, 
and smooth solid lines on the right-upper plots in Figures [2] and [31 The fastest growing 
unbound eigenmode grows at a rate Aq — 0.3323, and its spectral eigenfunctions are shown 
by the red jagged spiky lines. Next, in the case Pm ~ 6.7 x 10~^ and h = 1 there are 
eighteen bound magnetic eigenmodes in total, all are shown on the plot (C) in Figure [TJ 
The spectral eigenfunctions of four selected bound eigenmodes Ai ~ 41.73, A2 — 49.56, 
Aio — 695.9, and Aig ~ 17055 are shown by the dotted, dash-dotted, dashed, and smooth 
solid lines on the left-lower plots in Figures [2] and [31 The spectra of the fastest growing 
unbound mode Aq — 39.57 are again shown by the red jagged spiky lines. Finally, in the case 
Pm ~ 6.7 X 10~^ and h = 0.1 there are four bound eigenmodes, Ai ~ 0.39581, A2 — 12.30, 
A3 ~ 103.5, and A4 ~ 669.8, refer to the plot (D) in Figure [H These four modes and the 
fastest growing unbound mode Aq — 0.39573 have spectra that are shown on the right-lower 
plots in Figures [21 and [31 by the dotted, dash-dotted, dashed, smooth solid, and red jagged 
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spiky lines, respectively. 

Based on the results presented in Figures [TH31 we make the following important obser- 
vations. 

First, when the kinetic helicity increases, the number of bound magnetic eigenmodes 
increases significantly. Their growth rates, A„, become strongly concentrated near the growth 
rate of the fastest unbound eigenmode, Xq. (This last result follows from a nearly uniform 
distribution of A„ on the logarithmic-scale plots (A) and (C) in Figure [H) 

Second, on all plots in Figure[T]the growth rate of the first bound eigenmode, Ai, happens 
to be very close to Aq. We found same result for all other high Reynolds number cases that 
we investigated (not reported here), with different admissible values of the helicity parameter 
h. Thus, we propose that for dynamo driven by high Reynolds number Kolmogorov-type 
velocity field there always exists a shallow bound eigenmode Ai, such that Ai — Aq ^ Aq. 
This shallow mode grows faster than any of the unbound modes because Ai > Aq. 

Third, the spectra of the shallow bound eigenmode Ai, shown by the dotted lines in 
Figures [2] and [31 and the spectra of the fastest growing unbound eigenmode Aq, shown by the 
red jagged spiky lines, are close near the magnetic energy containing large scales. (Small- 
scale structures of these modes are however different: the shallow eigenmode has a relatively 
larger small-scales component.) In practical applications, the shallow mode grows faster 
than the unbound modes and can dominate at large scales; however it cannot be described 
by the conventional a-model for large-scale dynamo, since this model does not capture the 
bound modes. 

Fourth, consider the spectra of eigenmodes Aq and Ai (the dotted fines and the red 
jagged spiky solid fines in Figures [2] and [3]) . When the value of the kinetic helicity drops 
by a factor of 10 (from h = 1 to h = 0.1), the location of the peaks of these spectra shifts 
to larger scales by the same factor. Thus, the characteristic scales of eigenmodes Aq and 
Ai are both approximately equal to ~ lo/h, so that both these modes peak at large scale 
when kin etic helicity is small. % This resu lt is consistent with general predictions of the 



a- model (ISteenbeck. Krause fc Radlerlll966l ) 



At large correlation scales x 
Ao - .g2(0)/[KL(0) + 277] ^ h^vo/lo, 



^ 00, the eigenfunction of the fastest growing unbound eigenmode, 
asymptotically becomes a mixture of cosine and sine standing waves 
h/ln, whil e the bound eigenfunctions decline as oc exp(— / i„.a:), 



with wavenumber ko — \/An / a/k (0) + 2r7 
where fi„ = A„ — Aq / \/kl(0) + 27] (|Boldvrev. Cattaneo fc Rosner.,2005; .Malvshkin fc Boldyrevll2007n . 
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4. Conclusion 

We have presented the full numerical characterization of the kinematic Kazantsev- 
Kraichnan dynamo model in the case of the Kolmogorov scaling of the velocity field. Our 
main conclusion is that the structure and the characteristic correlation scales of the magnetic 
field amplified by helical kinematic dynamo action are determined by both bound (localized) 
and unbound (nonlocalized) growing eigenmodes. In particular, the large-scale component 
of the field is defined by the unbound eigenmodes and by the shallow bound eigenmodes. 
Because these shallow bound eigenmodes have growth rates higher than those of all unbound 
eigenmodes, at any given scale the former may rapidly become dominant over the latter. In 
practical applications, this means that the shallow bound modes, rather than the unbound 
modes, are likely to become essential in the large-scale magnetic field configurations in astro- 



physical systems. In this case the conventional a-dynamo model (ISteenbeck. Krause fc Radler 



19661 : lMoffattlll978l : lKulsrudll2005l ) gives an inadequate description of the large-scale magnetic 



field even at the kinematic stage of dynamo action. 

The a-model becomes inapplicable in this case because it uses a critical assumption that 
small-scale fiuctuations of the velocity and magnetic fields are much weaker and concentrated 
at the scales much smaller than the scales of the growing large-scale field. Under this 
assumption the a-model is obtained by averaging the induction Equation ([T]) over these 
small-scale fiuctuations to obtain a linear and homogeneous differential equation for the 
large-scale mean magnetic fieldB Thus, the a-model misses all growing bound magnetic 
eigenmodes, including the essential shallow bound eigenmodes that determine the eventual 
large-scale configuration of the magnetic field. 

We thank Fausto Cattaneo for many useful and stimulating discussions. This work was 
supported by the NSF Center for Magnetic Self- Organization in Laboratory and Astrophys- 
ical Plasmas at the Universities of Chicago and Wisconsin-Madison. S.B. is supported by 
the U.S. Department of Energy under the grant no. DE-FG02-07ER54932. 
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Fig. 2. — Absolute values of magnetic energy spectral eigenfunctions for four selected bound 
eigenmodes (shown by the dotted, dash-dotted, dashed, and smooth solid lines), and for 
the fastest growing unbound eigenmode (shown by the red jagged spiky solid lines). The 
left-upper, right-upper, left-lower, and right-lower plots are for the cases Pm ~ 150 & /i = 1, 
Pm ~ 150 & /i = 0.1, Pm ~ 6.7x 10"^ & /i = 1, and Pm ~ 6.7x 10"^ Lh = 0.1, respectively. 
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Fig. 3. — Same as Figure [2] except the absolute values of electric current helicity spectral 
eigenfunctions are plotted here. 



